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Abstract 

Atoms confined in a magnetic trap can escape by making spin-flip 
Majorana transitions due to a breakdown of the adiabatic approxi- 
mation. Several papers have studied this process for atoms with spin 
F = 1/2 or F = 1. The present paper calculates the escape rate for 
atoms with spin F > 1. This problem has new features because the 
perturbation AT which allows atoms to escape satisfies a selection rule 
AF Z = 0,±1,±2 and multi-step processes contribute in leading order. 
When the adiabatic approximation is satisfied the leading order terms 
can be summed to yield a simple expression for the escape rate. 

PACS numbers: 32.80.Pj, 41.20.Gz, 03.65.-w 

1 Introduction 

This is a more complete version of a Brief Report published in Physical Re- 
view A 74 035401 (2006). Magnetic traps for neutral atoms have many ap- 
plications in atomic physics and quantum optics. For example they are used 
to confine and study Bose-Einstein condensates. Magnetic wave guides can 
transmit and manipulate atomic de Broglie waves. Atoms may be trapped in 
fine structure level or hyperfine levels depending on the nuclear spin. Alkali 
atoms which have a non-zero nuclear spin are trapped in hyperfine states, 
while atoms of 52 Cr (Schmidt et al [1]), which have a zero nuclear spin but 
a large electronic spin are trapped in fine structure states. In this paper we 
focus on atoms trapped in hyperfine states, but the theory of fine structure 
trapping is essentially the same. 

Atoms are confined only in certain Zeeman levels and spin-flip transitions 
cause them to be lost from a trap. Thermal fluctuations in the surrounding 
apparatus can cause atoms to make spin-flip transitions to untrapped states. 



The theory of thermally induced losses has been developed by Henkel et al [2] 
and has been studied experimentally by Harber et al [3]. Losses due to radio- 
frequency noise in the currents that form the microtrap have been studied 
by Leanhardt et al [4]. Other references can be found in these papers. Spin- 
flip transitions between Zeeman levels due to a breakdown of the adiabatic 
approximation can also cause atoms to escape from a trap. These transitions 
were first studied in a one dimensional time-dependent model by Majorana 
[5] in (1932) and are often called Majorana transitions. In many experiments 
Majorana losses are much less important than those due to environmental 
effects; but they do come into play in some experimental situations (Schmidt 
et al [1], Ott et al. [6]). 

The spin of an atom in a magnetic trap precesses about the direction of 
the local magnetic field with frequency uj prec . At the same time its center of 
mass oscillates in the field with frequency uj V ib- In the adiabatic limit, when 
^vib « Uprea the component F z of hyperfine spin F along the direction 
of the local field is an approximate constant of the motion. Under these 
conditions the atom moves in an effective potential V^(r) = /j,ogF z B(r) 
where B(r) is the magnitude of the local magnetic field B(r), /j,q is the Bohr 
magneton and g is the hyperfine g-factor. Spin states with gF z > can be 
trapped near a minimum of B(r). 

Corrections to the adiabatic approximation [7] cause spin-flip a transition 
to a state with gF z < and the atom escapes from the trap. The cases F = 
1/2 or 1 were studied in Sukumar and Brink [7] and the transition rates for 
atoms to escape from a trap were calculated. Similar calculations have been 
made by other authors [8]. In the present paper we extend the results of 
ref.[7] to cases where the hyperfine spin F > 1 

The perturbation terms AT\ and AT2 in eq.(5), which cause a break- 
down of the adiabatic approximation, are calculated in section 4 and the 
transition rate for spin flip transitions is calculated in section 5. The dom- 
inant perturbations contribution of AT\ has a selection rule AF Z = 0, ±1 
and, for F > 1 while AT2 also allows transitions with AF Z = ±2. When 
F > 1 higher order perturbations and interference effects are important. 

Explicit expressions for Majorana transition rates are given in section 5. 
They show how the rates depend on the vibrational frequency of an atom 
in a magnetic trap, the Zemann splitting and the spin of the initial state. 
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2 The adiabatic Hamiltonian 



An atom with mass m, hyperfine spin F and magnetic moment \x moves in 
a magnetic field B(r). The interaction energy is 

V = —fx ■ B = nogF.H (1) 

where fio is the Bohr magneton and g is a ^-factor for the hyperfine state 
and the Hamiltonian for the motion of the atom coupled with the dynamics 
of the spin is 

2 

H = ^- + fi gF.B = T + h(r). (2) 
2m 

The adiabatic method starts as in ref.[7] by introducing basis states 
\F z (r)) which are eigenstates of the component of F along the direction of 
the local magnetic field. They are obtained from the eigenstates \F Z ) of F z 
by an r-dependent rotation 

|F,(r)> = U(r)\F z ). (3) 
The Hamiltonian in the adiabatic basis is 

H' = U- 1 (r)HU(r) = + V ad {r) + AT (4) 

2m 

where V a( i(r) is the adiabatic potential and AT is a correction coming from 
the transformation to the adiabatic basis. The adiabatic approximation 
corresponds to neglecting AT and corrections are obtained by treating it as 
a perturbation. The analysis of ref.[7] shows that AT = ATi + AT2 where 

ATi = i [A(r) - p + p - A(r)] ' AT2 = "L A - A (5) 

where A(r) = -^C/ _1 VC/(r). 

Following Sukumar and Brink [7] the unitary operator U(r) is chosen to 
be a rotation through it about an axis n(r) where n is a unit vector which 
bisects the angle between the z-axis and the direction of B(r) of the field 
when the atom is at r. Hence U(r) = exp(i7rn(r) • F). If the atom moves 
then n(r) —> n(r + Sr). To first order in Sn 

U^SU = 2in x Sn ■ F. (6) 

The resulting expression for A(r) is 

A(r) = -ihU\r)[VU(r)] = 2hn x Vn • F. (7) 
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In the adiabatic basis F has components F z along the direction of the 
local field B and the adiabatic potential is 

V ad (r) = liogF z B(r) (8) 

where B(r) is the magnitude of the magnetic field. There is a set of adiabatic 
potential surfaces, each with its own value of F z . For simplicity we assume 
that flog > so that states with F z > can be trapped around the minimum 
of B(r) at r = r m . When F z < the adiabatic potential V a d has a maximum 
at r = r m and the states are not trapped. 

Trapped atoms in states with F z > can escape from the trap by making 
a transition to a state with F z < because of non-adiabatic effects. The 
perturbation AT mixes states with different F z and can cause transitions 
out of the trap. The quantity A(r) is proportional to the spin F so that the 
perturbation AT\ can mix states with AF Z = 0, ±1. The perturbation AT2 
has a selection rule AF Z = 0, ±1, ±2. Trapped atoms with F z = or 1 
can escape by making a single step transition. This was the case discussed 
in ref. [7]. For atoms with F z > 1 higher order effects must be included. 



3 A simplified choice for the magnetic field 

The results in the last section hold for a general magnetic field. Many 
different devices for trapping neutral atoms have been developed. A typical 
configuration is the Ioffe-Pritchard trap [9]. The fields are produced by 
current carrying coils or by permanent magnets. Near the center of the trap 
the magnetic field is approximately 

B^B k + B' p (xi - yj) + ^ [z 2 - k + ... (9) 

and is characterized by an axial bias magnetic field Bq, an axial field curva- 
ture B z " and a radial quadrupole field gradient B p . 

Traps for ultracold atoms have been constructed on several different 
length scales and in many of these cases the field near the center of the 
trap can be approximated by eq.(9). The fields in conventional macro- 
traps, like the one used in the experiments of Harber et al [3], have a size 
on the scale of several centimeters. Typical values for the radial gradient 
and axial curvature are 50 G/cm and 100 G/cm 2 respectively. The fields 
in micro-traps are produced by current carrying wires on the surface of a 
microchip (Folman et al [10]). Distances from the wire to the center of the 
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trap are in the range 100-1000 fim with a radial field gradient in excess of of 
10,000 G/cm. Atom wave guides have a non-zero axial bias field and radial 
gradient but with an axial curvature B z = so there is no confinement in 
the z-direction. 

To simplify the calculations in this paper the field B(r) is assumed to 
have a constant axial bias field Bq and a radial field gradient B p = A. The 
axial curvature is set equal to zero so that there is no confinement in the 
z-direction. This can describe the field in an atomic wave guide with its 
axis parallel to the z-axis or a simplified trap where the atom is restricted 
to move in the (x, y) plane. The spin has components F z along the direction 
of the local field B and the adiabatic potential is 



Vad = vogF z B(r) = wgF^Bl + X 2 (x 2 + y*). (10) 

When the uniform field Bq is large enough the adiabatic potential can 
be approximated by a harmonic oscillator 

V ad « pogFz (b + ^-(x 2 + y 2 )^ . (11) 

The adiabatic potential depends on F z and both the oscillator frequency 
ijj v i\) and length parameter b depend on F z . They are given by 



u = uj^F z , b 2 = b\l^F z 



where 



mBo ' mujQ y m/iogX 2 

We assume that flog > so that states with F z > are trapped in the 
adiabatic potential. When F z < then Vad is an inverted parabola and the 
states are not trapped. 

The energy gap between the adiabatic surfaces at r = is 

Eq = HogBo = hiO prec (13) 

where ujp rec is the precession frequency of the spin in the magnetic field at 
r = 0. The adiabatic approximation is reliable when 



_ uj _ Hloq _ \ 2 bl _ I h 2 X 2 i 

X0 Up rec E Bq V m ^0g B 
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This is equivalent to the adiabatic condition in [7]. The parameters uq, Bo, 
Eq and xo depend only on the fields and not on the spin. 

The initial state of the atom in the trap has spin component F z = F z i 
oscillator length parameter 6, and oscillator frequency uii . The energy of 
the atom in the initial state is 

Ei = FzifiogBo + TiLQi ps F z iLo prec (15) 

when xo << 1- 

In the case of integer spin there the adiabatic potential is zero in the final 
state and the atom has a kinetic energy fixed by the energy conservation 
condition 

Ei = F Z iHogBa = TiLo prec F Z i. 

We see from the second last term in eq.(14) that the condition that 
the frequency ratio xo is small ensures that the expansion (11) is a good 
approximation for the low states in the adiabatic potential. 



4 Calculation of A7\ and AT 2 

The spinor field A(r) is given by eq.(7) where the unit vector n is chosen as 
in ref.([7]) 

n = (3xi — Pyj + ak (16) 



with 



2 Eq + B ^ 9 A 



" = -is-- " = wiB^wy <17) 



Substituting the expression for n x Vn • F into eq.(7) gives 

A(r) = Ai(r) + A 2 (r) + A 3 (r) with (18) 

Ai(r) = 2ha(3(F y Vx + F x Vy), (19) 

A 2 (r) = 2hf3 2 (yVx-xX7y)F z , (20) 

A 3 (r) = 2h(yF x + xF y )(aVP- pVa). (21) 

This equation for A(r) corresponds to one given in ref.( [7]). The explicit 
expression for (3 2 is given in eq.(17) and related factors in the equations for 
Ai(r) and A 3 (r) simplify to 

A A 3 r 

af3 = —, (aV/3-/3Va) = 



25' v ' 2B 2 (B + B )' 
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The eigenstates in the adiabatic potential are \L,n,F z ), where L is the 
orbital angular momentum about the symmetry axis and n is the number of 
nodes in the radial wavefunction. The energy eigenvalues in the harmonic 
approximation are 

E = n gF z + y/F^hu) (\L\ +n + l). (22) 

The matrix elements of AT conserve the quantity L — F z . This can be seen 
explicitly by examining the contributions of the various terms in ATi and 
AT2 and also, more generally, from the symmetry properties of the magnetic 
interaction which is proportional to 

B F = B F Z + xF x - yF y = B F Z + -{x + F + + x_F_) (23) 

where x± = x ± iy and F± = F x ± iF y . 

The terms A2 • p, Ai • Ai and A2 • A2 in AT all commute with F z while 
A2 • A3 = 0. These terms do not contribute to transitions out of the trap. 
The leading contributions in the adiabaticity parameter xo come from the 
terms proportional to Ai ■ p and A x • A3. They are 

Fl = i Al • p = ^B^ p y = -*4^+ F + - p - F -^ (24) 

and 

2h 2 A 4 h 2 A 4 

(25) 

where p± = p x ± ip y . The contribution of V\ in second order and the 
contribution of V2 in first order are both proportional to Xo- I n the following 
B is approximated by Bq. 



5 Transitions out of the trap 

Atoms with a g-factor g > and with F z > are trapped in the adiabatic 
potential when the perturbation terms ATi and AT2 are neglected. They 
can escape by making a transition to states with F z < under the influence 
of the perturbations. The dominant terms in a perturbation series are V\ oc 
Ai -p/2m with selection rule AT 2 = ±1 and V2 oc Ai • A3 /2m with selection 
rule AF Z = 0, ±2. Trapped states with F z = 1/2 or 1 can escape by a 
single step transition under the influence of V\ . This case was discussed in 
ref. [7]. 
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In the following section it is assumed that the initial spin state is the 
ground state in the adiabatic potential with F z = F z i, radial quantum num- 
ber n = and angular momentum component L z = L = 

\L,n,F z ) = \0,0,F zl ). 

The perturbation AT introduces admixtures of states with smaller values 
of F z and values of L satisfying the selection rule A(F Z — L) =0. The 
emission of the atom from the trap depends on the amplitudes of states with 
F z = 2, 3/2, 1, 1/2 in the initial state. These components allow transitions 
to unbound states with F z < under the influence of the perturbation and 
the initial state acquires a width T. When the adiabatic approximation is 
valid r < huji < fiLO prec where is the vibrational frequency in the initial 
adiabatic potential and u> prec is the precession frequency of the initial state. 

The orbital part of the initial state is an s-state harmonic oscillator wave 
function 

^ (r) = ^ exp (-w) 

in the adiabatic limit. The leading contributions to the transition rate out 
of the trap are due to the parts V~ and Vf of the operators V\ and V2 
which decrease F z . They are 

where 60 is defined in eq.(12) and xo m eq.(14). The perturbation expression 
for the transition amplitude has contributions from and as well as 
interference terms. 

The initial state with L = 0, n = and F z = F z i is denoted by \(j>i), the 
final continuum state with F z f = or —1/2 by |/) and states with L = —j 
and F z = F z i — j by \j, n). The number of positive values of F z is denoted 
by p = F z i — F z f . The amplitude A for the decay out of the trap is given by 
the perturbation formula (Messiah, Chapter XVI, section 6) 



{f\V \J m ,n m ) — 

'jm ,n m ^l- 



(j2,n 2 \V \ji,ni) (ji,ni\V |0) 

E2,n 2 — Ei E\, ni - Ei 



(26) 
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where V~ = V{~ + V 2 ~. Each matrix element in the series (26) contains 
either V{~ or V 2 ~ , but not both, because they satisfy different selection rules 
for AL. The general term in the series (26) contains p\ matrix elements of 
V{~ and P2 matrix elements of V 2 ~ where p\ + 2p 2 = p is the 'number of 
perturbation steps; 

p = F z i for integer spin 

P = Fzi + 1/2 for half integer spin. (27) 

In equation (26) Ei is the energy of the initial state and Ej mt n m is the 
energy of the state j m ,n m . They are given by 

Ei = E()F z i + huo \J F z i , Ejm,n m = Eo(F z i—j) + (2n m + l)huo \J F z i — j m 

and all depend on F z because the adiabatic potential depends on F z . The 
expression for A is complicated because there is mixing of many orbital 
states rii for each ji. 

There is a big simplification in the adiabatic limit. Then Tito « Eq and 
the energy denominators can be approximated by 

E jm,n m - -Ei ~ -jEo 

and the sums over m, ri2, ■■■ in the perturbation formula for A\ can be eval- 
uated by closure. The structure of the final expression for the transition 
amplitude is 

A^fu^oh^X XQ ^(F zf \Fl\F zi ) ]T W-1) P2 AW> 2 (28) 

pi+2p 2 =P 

where 

W = ^yih^-^i) = [y) (-lY»{kf,L f \(bop-/Ky\<lH) 

(29) 

and kf is the wave number of the final free article state. The phase factor 
(— 1) P2 in eq.(28) arises because all the energy denominators in (26) are 
negative, and the one in eq.(29) comes from the reduction of the matrix 
element. They cancel and all the interference terms between and V 2 ~ 
are constructive. The N PjP2 are numerical combinatorial factors. For p2 = 
and p2 = 1 they are given by 

Some values for p2 = 1 are given in the Appendix. 
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5.1 Integer spin 

The expression for the transition rate out of the trap depends on whether 
the hyperfine spin of the atom is integral or half integral. In the case of 
integral spin the adiabatic potential is zero in the final state and the orbital 
state of the atom is a plane wave 4>f = {L)~ l exp(iky • r) normalised in 
a box with side L. The magnitude of ky is fixed by energy conservation 
k"j = 2mEi/h 2 . The orbital matrix elements can be evaluated and 



where 



1 Z,l 



io(k/) = — e t l/ (30) 
is the overlap of the initial and final orbital states. Hence 

A « huo {^^\t l {F zf \F p _\F zl )^(k x - ik y ) p I (k)C p . (31) 



The factor 




N P , P2 ■ (32) 



Some values of N PjP2 / {F z ,j) P2 and C p are given in Appendix A. The numbers 
show that the contribution of the perturbation AT2 = A.A/2m to the 
transition rate is important. This is because the terms in the perturbation 
formuala add coherently. As explained in section 5 p is related to the spin 
component of the initial trapped state by p = F z i for integer spin and 
P = F z i + 1/2 for half- integer spin. 

The transition rate out of the trap is given by Fermi's Golden Rule 

w f = —\A\ Pf where p f = -^-^ (33) 

is the density of final states. The square of the magnitude of the transition 
matrix element reduces to 

\A? « S^f {^fj" m ,\ F ^cf-^e-W (34) 

where we have used 

ul = uf/F zi and x^\k) = 2F zi . (35) 
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The first of the expressions in eq.(35) relates the frequency uji of the initial 
state to the standard frequency ujq defined in eq.(12) and the second comes 
from the energy conservation condition 

h 2 1 h 2 

—k) = F zi E = F zi huj — = F zi —2 — . (36) 
2m xo mboXo 

in combination with eqs.(12) and (14). Substituting into eq.(33) gives an 
expression for the escape rate from the trap 

w * = lr (¥)" 1 i^/^-i^)! 2 ^ 2 ^- ( 37 ) 

where p = F z i (eq. 27). 

In the case of integral spin the spin of the atom in the unbound final 
state is F z f = and the angular momentum factor is 



(38) 



|(0|^|^)| 2 = F(F + l)(F-l)F + 2)..(F-p + l)(F + p) = ||±^: 

where F is the total spin of the atom. 

The exponent in eq.(37) can be written in a number of different ways by 
using the relations 

k 2 b 2 = hllk = 2 -^. (39) 



Equation (37) shows that, for a given kf, the transition rate is reduced 
when the width 6j of the initial state decreases. The second term in eq.(39), 
(kfbi) 2 = 2y/F z iXo shows that the escape rate is exponentially small when 
the adiabaticity parameter is small (xo << 1)- For a given xo it is reduced 
by an increase in the spin component F z i of the initial state. The pre- 
exponential factor Xo gives a further reduction in the escape rate but the 
factors depending on F and F z i increase the escape rate. Using {kfbi) 2 = 
2Ei/huJi the escape rate is proportional to exp(—2Ei/hu>i). It is small when 
the magnetic energy is large compared with the vibrational energy in the 
initial state (Ei » fiuji). 

Experimental papers normally give the axial bias field Bq in Gauss (G) 
and the radial field gradient B' = A in G/ cm. Expressions for uj prec and ujq 
in these units are 



uj prec = 8.8 x lQpgBo sec _1 G _1 , oj = 74.6A A /-^- sec _1 G _1 . (40) 

ABn 
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Values of the parameters vary a lot from one experiment to another. For 
example a microtrap described by Ott et al [6] for 87 Rb atoms operated in 
several different modes. In one mode -Bo ~ 0.7 G and A ~ 122 G cm -1 so 
that uj ~ 5x 10 3 sec _1 and u prec w 1.7 x 10 6 . The ratio xo ~ 2.9x 1CT 3 << 1 
and the adiabatic approximation is expected to be reliable. This situation 
is tyical of most Ioffe-Pritchard traps. 

There are cases where the adiabatic parameter xo ~ 1- This happens 
if the radial field gradient is very large (cf [6]) or if the axial bias field is 
very small (cf [1]. In the second mode of the trap described by Ott et al 
[6] B Rs 1 G and A « 3 x 10 5 G cm' 1 . In this case u Q « 3.7 x 10 6 sec" 1 , 
ujp rec ~ 2.4 x 10 6 and xo ~ 1-5 and the adiabatic approximation is not valid. 

The expression (37) is interesting because the exponential factor is pro- 
portional to probability density Po(kf) for finding a final momentum kf 
when the initial orbital state of the atom is the ground state in the adia- 
batic potential. The expression can be generalized to the case of a general 
initial orbital state with momentum distribution P(k); 

"f = ^ (^) P 1 \(F Z f\F P -\F zi )\ 2 C^P(k f ). (41) 

Here the C p are different from the values in Table 1 but are still of the 
order of unity. In this form the spin-flip transition rate is small when the 
probability for finding the momentum hkf in the initial state is small. 

5.2 Half integer spin 

When the hyperfine spin is F = 1/2, 3/2... the final value of the spin projec- 
tion is F z j = — 1/2 and the final orbital wave function is a continuum state 
in an inverted parabolic adiabatic potential. The changes in the expression 
for the transition rate can be estimated by replacing the plane wave final 
state by a semiclassical approximation 

Mr) « (i + r 2/( W/) ex P (/ \A 2 + ^/) 

which reduces to the plane wave near r = 0. Here kf is the momentum of 
the atom in the final state at r = and bf is the oscillator length parameter 
of the inverted parabolic potential. If the integrals in the matrix elements 
are estimated by the method of stationary phase then the expression (37) 
for the transition rate is replaced by 
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wf = im fpgY 1 \ { - 1/2 \F*_\F Zi )fcy ck2 f b l (42) 



2 

where p = F zi + 1/2 (eq.27) and 



c = V2F~t 3 ,n- 1 ( ] j^-). (43) 

The biggest effect of the factor c is for a transition from an initial state 
with F z i = 1/2 to a final state with F z f = — 1/2 when c = -it/ 4. For large 
F z i the factor c — > 1. (There was a printing error in eq.(34) of ref.[7]. A 
factor 7r was omitted in the exponent). 

The angular momentum factor is replaced by 

(-l/2|F5|F^ = (F + l/2)^±^|. 
With these substitutions eq.(37) reduces to 

ITLUi ( TTkffi\ 



w f = — exp 



when F = F zi = 1/2. 



6 Discussion and conclusions 

We have obtained approximate expressions for the decay rate of an atom 
from a magnetic trap due to breakdown of the adiabatic approximation. 
The results hold for an arbitrary initial spin state of the atom. When the 
initial orbital state of the atom is the ground state in the adiabatic potential 
the decay rate in eq.(37) is proportional to the vibration frequency of the 
initial state multiplied by an exponential factor (eqs. 14, 37, 39) 

exp(-dfejb?) = exp (-2c^p) = exp (^-2c^^j . (44) 

The factor c = 1 for integer spin and is given by eq.(43 for half-interger 
spin. There is also a factor (pXo/^) P_1 ' w here p = F z i for integer spin and 
p = F z i + 1/2 for half-integer spin. This further reduces the transition rate 
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in the adiabatic limit. These expressions generalize the results of [7] where 
p = l. 

Under normal experimental situations Majorana transition rates are 
small because the adiabaticity parameter xo is small and the ratio of spin 
precession frequency is similar to the basic trap frequency E /fr ~ uj is 
large. Under these conditions the exponential factor (44) is very small. 
Eq.(14)shows that xo °^ (Bo) ~ 3 / 2 . Thus Majorana transitions can be seen 
when the bias field B$ is reduced sufficiently to make xo ~ 1- The cal- 
culations in this paper use the harmonic approximation for the adiabatic 
potential. The arguments at the end of section 3 show that the harmonic 
approximation is valid whenever the adiabtic approximation is valid. 

After the escape of the atom from the trap it has a kinetic energy equal 
to Ei. The transition rate is small when the final kinetic energy is much 
greater than the zero point center of mass energy in the initial state or 
when kfbi >> 1. The argument leading to eq.(41) shows that the transition 
rate is proportional to the probability density Po(kf) for finding the final 
momentum kf in the bound state in the initial adiabatic potential. It be- 
comes small when kf is large. This form shows that decay rate out of the 
trap is largely determined by momentum matching between the initial and 
final state. At a high temperature Po(kf) is replaced by a Boltzmann factor 
exp(— h 2 k 2 f /2mk B T) as in ref. [7]. 

The analysis in section 4 shows that the perturbation terms proportional 
to A(r).p and A(r).A(r) both contribute in leading order to decay out of 
the trap. Their contributions interfere constructively when the orbital state 
is the ground state in the initial adiabatic potential. The interference effects 
are contained in the factor C 2 (eq. 32). 

The vibrational frequencies uj\\ and u)± of states in a 3-dimensional trap 
are determined by the axial field curvature B" z and the radial curvature 
A = B' p . A detailed discussion of losses from a 3-dimensional trap due to 
a breakdown of the adiabatic approximation is more complicated than the 
case considered here; but qualitative considerations based on eqs.(37) and 
(41) suggest that losses due to the parallel motion will be negligable when 

(J 1 1 << L0±. 
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Appendix: The coefficients N, 



P,P2 



This comes form the energy denominators in the perturbation formula. 
There are p — 1 energy denominators and N PtP2 = l/(p— 1)! when p2 = 0. 
When p2 = 1 one of the energy denominators is missing. This can happen 
in p — 1 ways. 

Simple counting gives 



1 



p,0 



N, 



(P-1)! 

p(p 



when P2 = 



1) 



2(p-l)! 



when P2 = 1 



The minus signs come because the energy denominators are all negative. In 
general it is (— 1) P2 . We could not find expressions for general p and P2 > 1. 
Numerical values for £>2 = 2 and p < 5 are given in the following table. 





p = 1 


p = 2 


p = 3 


p = 4 


p = 5 


N P,0 


1 


1 


1/2 


1/6 


1/24 


N P ,i/p 





1/2 


1/2 


1/4 


1/12 


N p , 2 /P 2 











1/32 


1/40 


C p 


1 


3/2 


1 


43/96 


3/20 



Table 1, Values of AT p o, -/V P ,i, AT P)2 and C p . 
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